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Preface 

On  two  occasions,  the  U.S.  Air  Force  research  and  development  community  has  proposed  multi¬ 
function  space-based  sensor  concepts  that  employ  a  formation  of  identical,  free-floating,  semi- 
autonomous,  and  self-organizing  micro-satellites  to  form  large  space  apertures  [1,2].  In  view  of 
recent  advances  in  space-time  processing,  we  consider  here  the  design  principles,  performance, 
and  technical  risks  associated  with  such  a  network  of  satellites  that  function  collaboratively  as  a 
distributed  array  radar  for  detecting  and  tracking  moving  objects  near  the  earth  surface. 

The  satellite  network  is  viewed  as  a  cost-effective  approach  to  achieve  the  large  apertures  in  space 
required  for  detecting  slowly-moving  targets  with  a  rapidly-moving  satellite  array.  Cost 
effectiveness  is  envisioned  from  efficient  mass  production  of  nearly  identical  micro  satellites, 
opportunistic  launch  and  gradual  network  formation,  and  maintenance  and  performance 
enhancement  through  gradual  replacement  of  network  satellites  during  mission  lifetime  with 
residual  space  launch  capacity.  In  addition,  the  satellite  network  can  conceivably  accomplish 
multiple  missions  via  reconfiguration  to  provide  passive,  interferometric  emitter  location,  local 
and  wide  area  communications,  and  a  variety  of  radar  functions  including  synthetic  aperture  radar 
mapping  and  detection  and  tracking  of  targets  moving  with  respect  to  the  earth  surface.  In  this 
work,  we  analyze  the  wide-area  radar  detection  of  slowly  moving  targets,  perhaps  the  most 
technically  demanding  of  the  potential  radar  missions  for  the  distributed  sensor  array. 

There  are  myriad  technical  problems  associated  with  distributed  networks  of  autonomous 
satellites  that  must  function  collaboratively  to  accomplish  a  common  mission.  Formation  flying  is 
one  problem.  The  satellites  must  navigate  in  the  formation  with  knowledge  of  the  common  array 
mission  and  information  on  the  locations  of  other  satellites  in  the  formation  to  avoid  collisions 
and,  desirably,  to  obtain  proper  formation  positions  with  adequate  precision  during  the  mission 
observation  time.  Navigation  must  be  achieved  with  minimum  energy  consumption  to  extend  the 
life  of  network  satellites.  A  second  problem  is  achieving  collaborative  behavior  of  the  satellites 
toward  the  common  array  mission.  The  architecture  for  centralized  mission  planning  with 
distributed  satellite  control  must  be  refined.  High  speed  inter-satellite  data  communications  are 
necessary  to  ( 1 )  effectively  locate  each  satellite  in  space  from  inter-satellite  range  measurements, 
(2)  achieve  common  transmit  and  receive  signal  timing  and  coherence  for  sparse  array 
beamforming  and  signal  processing,  and  (3)  implement  transmit  waveform  selection  and 
scheduling  and  array  processing  for  beamforming  and  receive  signal  processing.  Finally,  a  third 
problem  is  associated  with  the  performance  of  a  search  radar  operating  against  ground  reflections 
or  clutter  when  implemented  with  a  large,  coherent,  highly-thinned  array  on  transmit  and  receive. 
The  large,  sparse,  multi-element  time-varying  phased  array  can  form  a  narrow  beamwidth  array 
pattern  with  concomitant  grating  and  random  sidelobes  that  introduce  significant  ground  clutter 
into  the  receive  signal. 

Previous  work  described  approaches  to  distributed  space-based  radar  employing  satellites  in  a 
single  orbit  [3-7].  This  report  summarizes  our  work  [8-1 1]  on  a  more  advanced  system,  which 
employs  a  cluster  of  free-floating  satellites  in  a  unique  periodic  constellation.  The  geometry  is 
such  that  orbital  mechanics  maintains  cluster  periodicity  as  the  cluster  orbits  the  Earth.  The 
satellites  have  identical  phased  array  antennas  for  transmit  and  receive,  and  each  satellite 
transmits  its  own,  possibly  orthogonal,  signal  and  receives  the  reflected  signals  from  all 
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transmitting  satellites,  similar  to  the  French  RIAS  groundbased  system  [12-14].Thus  the  system 
constitutes  a  full  multiple-input  multiple-output  (MIMO)  system. 

Later  we  investigated  signal  processing  techniques  and  GMTI  performance  of  more  general  linear 
and  planar  random  arrays  [15-17].  These  results  will  be  published  in  a  companion  report  [18].  For 
a  linear  array  that  moves  along  its  axis,  it  is  shown  that  the  optimum  space-time  processing  leads 
to  a  solution  that  is  similar  to  a  displaced-phase-center  antenna  (DPCA). 
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Summary 

The  TechSat2I  space-based  radar  concept,  suggested  by  AFRL/VS,  employs  a  cluster  of  free- 
floating  satellites,  each  of  which  transmits  its  own  orthogonal  signal  and  receives  all  reflected 
signals.  The  satellites  operate  coherently  at  X-band.  The  cluster  forms  essentially  a  multi-element 
interferometer  with  a  concomitant  large  number  of  grating  lobes  and  significant  ground  clutter. 
This  study  explores  some  basic  characteristics  of  such  a  system  with  respect  to  detection  of 
slowly-moving  ground  targets  (GMTI). 

A  companion  report  [18]  investigates  the  GMTI  performance  of  basic  random  sparse  arrays.  For 
a  linear  array  which  moves  along  its  axis,  it  is  shown  that  the  optimum  space-time  processing 
leads  to  a  solution  that  is  similar  to  a  displaced-phase-center  antenna  (DPCA). 

In  this  report  we  propose  a  novel  two-dimensional  periodic  array  geometry  for  TechSat21  and  an 
associated  technique  for  pattern  synthesis  in  angle-frequency  space.  It  exploits  the  double 
periodicities  of  the  grating  lobes  in  the  angular  domain  and  of  the  radar  pulses  in  the  frequency 
domain,  and  allows  substantial  gains  in  clutter  suppression. 

We  consider  a  ‘separable’  system  architecture  for  TechSat21,  where  a  planar  M-element  array  is 
followed  by  a  single  N-pulse  Doppler  filter,  and  evaluate  some  processing  alternatives  for  their 
effectiveness  in  clutter  suppression.  An  attractive  feature  of  this  architecture  is  that  it  allows 
simple,  multi-pulse  MTI  processing  and  adaptive  antenna  pattern  control.  The  limited  number  of 
degrees  of  freedom,  M+N,  requires  a  relatively  low  computational  load.  General  space-time 
processing  is  also  evaluated  and  shown  to  offer  substantially  improved  signal/clutter  ratios. 
However,  the  larger  number  of  degrees  of  freedom,  MxN,  implies  a  heavy  computational  load. 

Finally,  taking  the  ideal  periodic  array  as  a  reference  case,  we  analyze  the  performance 
degradation  with  array  pitch  (satellite  cluster  rotation)  and  yaw  (Earth  rotation),  and  known 
random  element  position  errors.  It  is  found  that  the  signal/interference  improvement  factor  is  an 
extremely  sensitive  function  of  the  pulse  repetition  frequency  and  that  it  varies  over  a  large 
dynamic  range.  Achievable  improvement  factor  maxima  are  moderately  degraded  by  array  pitch 
and  strongly  degraded  with  increasing  position  errors.  Yaw  effects  appear  to  be  negligible.  In 
general,  with  increased  randomness,  fewer  elements  contribute  to  form  DPCA  pairs  so  that 
satellite  usage  efficiency  (and  SNR)  is  degraded. 

Thus,  a  periodic  array  lattice  allows  for  significantly  better  clutter  suppression  than  a  non-periodic 
lattice,  independent  of  the  particular  space-time  processing  scheme. 
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1.  Introduction  of  the  TechSat21  Space-Based  Radar  System 

The  space-based  radar  concept  TechSat21  [2,19,20]  is  based  on  a  distributed  cluster  of  free- 
floating,  small  satellites  operating  cooperatively  to  perform  a  surveillance  mission.  The  primary 
focus  is  the  Moving  Target  Indication  (MTI)  radar  mode.  A  baseline  assumption  is,  similar  to  the 
French  RIAS  system  [12-14],  that  the  individual  satellites  transmit  different  frequencies  (or  other 
orthogonal  signals),  but  receive  all  the  reflected  signals  at  the  frequencies  transmitted  by  the 
other  satellites  in  the  cluster,  as  shown  in  Fig.  1.  In  addition,  it  is  assumed  that  the  positional 
accuracy  can  be  sensed  to  within  a  fraction  of  a  wavelength  so  that  signals  can  be  combined 
coherently.  The  apparent  loss  of  array  gain  on  transmit  does  not  degrade  radar  target  detection,  as 
explained  in  Appendix  A.  Each  satellite  has  roughly  a  4  m2  phased  array  operating  at  about  10 
GHz,  other  typical  parameters  are  4  to  20  satellites  per  cluster,  cluster  diameters  100  to  1000 
meters,  and  orbits  700  to  1000  km  above  Earth. 

Our  companion  report  [  1 8]  considers  moving  target  detection  with  a  linear,  random  array  to 
derive  the  fundamental  characteristics  and  provide  insight.  In  the  present  report  we  summarize  our 
earlier  work  on  the  performance  of  a  more  specific  space-based  radar  system  with  a  sparse  two- 
dimensional  array  in  a  circular  low  Earth  orbit.  Other  approaches  to  distributed  space-based  radar 
employing  satellites  in  a  single  orbital  plane  are  described  in  [3,4,6,7],  We  discuss  a  particular 
array  geometry,  analyze  its  pattern  foot  print  on  a  spherical  Earth,  present  results  of  clutter 
suppression  for  signal  processing  techniques  of  various  levels  of  complexity,  analyze  the 
sensitivity  to  pulse  repetition  frequency  (prf),  and,  finally,  determine  the  effects  of  various  array 
errors. 


Figure  1.  The  TcchSat21  space-based  radar  employs  a  cluster  of  small  free-floating  satellites.  Each  satellite 
transmits  its  own  unique  signal  and  receives  all  reflected  signals  from  the  target. 


2.  Basic  Considerations 

Detection  of  a  slowly-moving  ground  target  by  a  moving  radar  is  based  on  the  difference  of  the 
Doppler  frequencies  of  the  target  return  and  the  ground  clutter.  The  situation  is  depicted  in  Fig.  2, 
where  the  ground  clutter  ridge /ciimCr  is  the  return  from  one  range  bin  illuminated  by  the  antenna 
main  beam  and  sidelobes.  The  different  angular  parts  of  the  range  bin  generate  an  extended 
Doppler  spectrum,  since  they  have  radial  velocities  that  vary  with  azimuth  angle  cp  as  seen  from 
the  radar.  A  moving  target  in  the  antenna  main  beam  appears  as  a  spike  in  the  angle-Doppler 
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domain.  To  separate  the  target  return  from  the  clutter  ridge  clearly  requires  a  narrow  antenna  main 
beam  and  a  narrow  Doppler  filter  bandwidth,  or  equivalently,  a  minimum  aperture  dimension  A 
and  coherent  integration  time  T.  These  are  shown  in  Appendix  B  to  be 

T  «  — — - —  ,  (1) 

^  ^ large/  ^  ^fargef 

where  V target  is  the  target  speed,  V  radar  the  radar  along  track  speed,  and  X  the  radar  wavelength. 
(Effects  of  the  illumination  angle  and  the  velocity  directions  can  increase  these  minimum 
requirements  but  have  been  suppressed,  as  discussed  in  Appendix  B).  Sampling  the  aperture  A  at 
A/2  spacing  leads  to  a  total  number  of  spatial  samples  (array  elements) 

/V  =3-^sl.  (2a) 

W 

target 

Similarly,  over  the  time  T,  the  total  number  of  temporal  samples 

N  -3^2L,  (2b) 

V 

v  target 

since  the  minimum  temporal  sampling  rate  equals  twice  the  highest  Doppler  frequency  2  VmdarlX. 

It  is  interesting  to  note  that  the  total  numbers  of  spatial  and  temporal  samples  are  equal,  and  that 
they  are  determined  solely  by  the  ratio  of  vmda/  v target,  independent  of  the  radar  frequency. 
However,  for  slowly-moving  targets  viewed  from  space  these  numbers  are  large.  For  example,  a 
target  speed  y target  =  2.2  m/s  (8  km/hour)  and  a  low  Earth  orbit  radar  speed  M  radar  =  7500  m/s 
leads  to 


Na=Nt=10\  (3) 

corresponding  to  a  (linear)  array  with  10000  elements  that  generates  108  space-time  samples  at 
each  range  sample  for  unambiguous  azimuth-Doppler  processing.  With  100  range  bins  and  a  prf  = 
1000  this  represents  a  data  rate  of  1013  samples  per  second.  It  is  obvious  that  the  narrow 
beamwidth  and  long  processing  times  can  slow  the  search  rate  of  such  a  radar  system  and  the 
large  number  of  space-time  samples  present  a  substantial  signal  processing  requirement.  Thus 
there  is  a  strong  motivation  for  reducing  the  number  of  samples,  i.e.  for  spatial  and  temporal 
undersampling. 
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Figure  2.  Target  return  and  ground  clutter  return  in  Angle-Doppler  space. 


A  spatially  undersampled  or  thinned  array  maintains  the  narrow  beamwidth  necessary  for  slow 
moving  target  detection  but  introduces  excess  clutter  associated  with  the  lost  sidelobe  control. 
Basically,  thinning  can  be  periodic  or  random,  which  leads  to  very  different  sidelobe  structures 
and  received  clutter  Doppler  spectra,  as  seen  in  Fig.  3.  Here  equal  transmit  and  receive  patterns 
are  assumed.  As  reference,  the  clutter  spectrum  corresponding  to  a  fully  filled  array  is  shown.  It 
decreases  monotonically  away  from  the  main  beam,  corresponding  to  the  continuous  decrease  of 
the  sidelobe  envelope.  Further  shown  is  the  spectrum  corresponding  to  a  periodically  thinned 
array,  whose  grating  lobes  generate  a  periodic  spectrum.  Here  the  clutter  is  concentrated  to 
relatively  high,  narrow  spectral  bands  with  low  clutter  regions  between.  Finally,  there  is  shown 
the  ensemble  average  clutter  spectrum  for  a  randomly  thinned  array.  The  clutter  power  is 
concentrated  at  a  Doppler  frequency  corresponding  to  the  main  beam  pointing  direction  with  a 
bandwidth  that  is  comparable  to  that  produced  by  the  fully  filled  array.  Outside  the  main  beam, 
the  power  is  uniform  at  a  level  proportional  to  1/N2,  where  N  is  the  number  of  elements  in  the 
array.  A  single  realization  of  the  randomly  thinned  array  is  also  shown  by  the  dotted  curve. 


RELATIVE  dbPPLER  FREQUENCY 

Figure  3.  Clutter  Doppler  spectrum  received  with  fully  Filled  array  (blue/green),  with  thinned  periodic  (blue) 
and  thinned  random  arrays  (red  solid  ensemble  average,  red  dashed  one  realization).  Two  way  propagation  is 
assumed. 


The  periodically  thinned  array  is  preferred  since  between  the  harmonic  spectral  lines,  the  power 
spectrum  of  the  clutter  is  lower  than  for  the  random  array,  permitting  better  target  detectability.  In 
addition,  it  is  compatible  with  temporal  periodic  undersampling  of  this  power  spectrum.  Selecting 
a  sampling  frequency  (pulse  repetition  frequency)  related  to  the  harmonic  spectral  lines  results  in 
aliases  which  overlay  the  original  spectral  lines,  and  thus  the  total  spectrum  retains  the  relatively 
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low  clutter  power  between  the  spectral  lines.  Fig.  4  illustrates  this  point.  It  shows  the  power 
spectrum  of  a  linear,  periodically-thinned  array  along  with  the  power  transfer  function  of  a  three- 
pulse  binomial  canceller,  and  the  resultant  uniform  residual  clutter  at  the  output  of  the 
cancellation  filter. 


Figure  4.  Incident  clutter  Doppler  spectrum,  Doppler  filter  transfer  function  and  resultant  output  power 
spectrum. 

The  array  signal  processing  considered  here  involves  an  array  beam  former  followed  by  a  Doppler 
filter  and  is  denoted  as  ‘separable’  processing,  since  it  performs  spatial  and  temporal  processing 
independently.  Fig.  5a  illustrates  this  separable  processing.  It  is  considerably  simpler  than  general 
space-time  processing,  where  each  space-time  sample  is  individually  weighted  and  summed,  (see 
Fig.  5b). 

This  basic  approach  of  a  thinned  periodic  array  with  a  ‘tuned’  pulse  repetition  frequency  (prf,  alt. 
/rep)  and  separable  processing  was  initially  evaluated  for  the  TechSat21  radar  system  [8-11].  In  a 
second  phase  we  evaluated  also  the  performance  with  full  space-time  processing  [15-16]. 


3.  Application  to  TechSat21 

Cluster  configurations  are  governed  by  orbital  mechanics,  i.e.  Kepler's  equations.  Linearization 
for  small  motions  around  a  reference  point  in  a  global  circular  orbit  leads  to  the  Hill  equations, 
derived  in  1878,  and  alternatively  named  Clohessy-Wiltshire  equations,  who  rederived  them  in 
1960  [20].  These  equations  constrain  the  satellites  to  move  either  along  a  linear  track  or  to  local 
orbits  around  the  reference  point,  such  that  their  projections  on  a  vertical  plane  form  2: 1  ellipses. 

Our  proposed  TechSat21  configuration  consists  of  a  vertical,  planar,  19-element  array  as  shown  in 
Fig.  6.  The  vertical  orientation  allows  looking  toward  both  sides  equally  well  and  gives  the 
highest  gain  at  maximum  range.  The  array  configuration,  with  one  central  element  and  six 
elements  on  each  of  three  concentric  2: 1  ellipses,  realizes  a  periodic  triangular  grid.  During  each 
Earth  orbit  the  elements  rotate  one  full  cycle  along  their  respective  ellipses  but  maintain  a 
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Figure  5.  (a)  Separable  signal  processing  -  array  beamformer  followed  by  Doppler  filter  provides  M+N  degrees 
of  freedom  for  Signal/Clutter  optimization,  (b)  Full  space-time  processing  provides  MxN  degrees  of  freedom. 


scan  -30°  -  60°  off  broadside 


Array  geometry 

19-ELEMENT  ARRAY  IN  TRIANGULAR  GRIO 


(a)  (b) 

Figure  6.  (a)  Vertical  array  in  Earth  orbit,  viewed  end-on.  Array  scans  30°-  60°  off  broadside,  (b)  Array 
Geometry:  19-Element  Array  in  Triangular  Grid. 


triangular  grid  (the  ellipse  major  axes  remains  tangential  to  the  circular  orbit).  The  triangular 
lattice  is  the  key  to  our  approach  since  it  generates  grating  lobes  that  are  highly  periodic  in  angle, 
and,  from  further  analysis,  periodic  in  Doppler  spectrum  as  well. 
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Although  simple  in  principle,  the  approach  is  complicated  by  the  curvature  and  rotation  of  the 
Earth.  The  Doppler  frequency  of  the  ground  clutter  is  derived  in  [18]  in  terms  of  the  satellite  orbit 
coordinates  and  the  array  look  direction,  and  the  result  is  repeated  in  Appendix  B  for  reference. 
Since  the  element  pattern  (the  pattern  of  an  individual  satellite)  illuminates  only  a  relatively  small 
angular  sector  of  the  Earth,  the  grating  lobes  and  Doppler  frequencies  all  experience  the  same 
linear  shift,  as  shown  by  numerical  analyses.  A  second  complication  is  the  rotation  of  the  array 
around  its  geometric  center,  which  continuously  changes  the  grating  lobe  structure.  However,  this 
can  be  accommodated  by  adjusting  the  pulse  repetition  frequency  accordingly. 


To  evaluate  the  above  approach  quantitatively,  we  developed  a  MATLAB  code,  which  computes 
the  received  ground  clutter  for  an  array  in  Low  Earth  Orbit.  Representative  results  are  shown  in 
Figs.  8,  9,  and  10,  which  apply  for  satellite  altitude  of  850  km  above  a  rotating  Earth,  orbital  plane 
inclined  70°  relative  to  the  equator,  and  a  satellite  position  45°  above  the  equator,  as  illustrated  in 
Fig.  7.  The  array  operates  at  10  GHz.  The  array  shape  is  elliptical  with  a  100m  major  axis,  and  a 
50  m  minor  axis.  The  19  elements  are  2x2  m~  phased  arrays,  the  array  rotation  angle  is  0°  and  the 
look  angle  is  45°  straight  down.  The  clutter  density  is  assumed  uniform  in  azimuth. 


Figure  7.  Illustration  of  the  orbital  geometry  chosen  for  the  numerical  examples 


Fig.  8  shows  the  mapping  of  three  range  bins  onto  the  array  pattern.  Due  to  its  particular  geometry 
the  range  bin  maps  into  a  straight  line,  u=const.,  in  the  array  pattern  when  using  the  so  called  sine- 
space  (u,v)-coordinates.  The  horizontal  center  line  represents  the  desired  range  bin  and  the  upper 
and  lower  lines  represent  aliased  range  bins,  corresponding  to  a  pulse  time  of  2  psec  and  pulse 
repetition  frequency  of  4640  Hz.  Note  that  the  aliased  bins  are  illuminated  with  very  low  intensity 
since  they  fall  outside  the  element  pattern  (one  satellite)  and  therefore  can  be  ignored.  The  desired 
range  bin  is  so  narrow  in  down-range  that  we  consider  the  pattern  constant  in  this  direction,  and 
only  consider  the  cross-range  variation  when  computing  the  power  of  the  Doppler  spectrum 
received  from  the  ground  clutter. 
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ARRAY  CONFIGURATION 


RANGE  BINS  MAPPED  ONTO  GRATING  LOBE  DIAGRAM.  prf=4640.  l-puls*=2*006 
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Figure  8a.  Configuration  of  19-element  sparse  array  and  contour  plot  of  the  array  pattern  (contours  at  -3dB,  - 
lOdB,  -20dB).  Horizontal  centerline  represents  desired  range  bin  mapped  onto  the  pattern,  upper  and  lower 
lines  represent  aliased  range  bins.  Array  rotation  angle  =  0  . 
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Figure  8b.  As  above,  but  array  rotated  through  an  angle  =  15°. 


Fig.  9  shows  the  main  part  of  the  received  Doppler  spectrum  corresponding  to  the  center  range 
bin  in  Fig.  8.  The  envelope  is  determined  by  the  two-way  element  pattern  (one  satellite)  and  is 
plotted  here  over  an  angular  sector  corresponding  to  the  main  beam  and  ±  3  sidelobes.  The  narrow 
‘spikes’  are  the  periodic  grating  lobes;  they  are  separated  by  approximately  1 160  Hz.  A  periodic 
pulse  train  illumination  creates  aliasing,  and  choosing  the  pulse  repetition  frequency  frepsuch  that 
the  4lh  grating  lobe  aliases  onto  the  main  beam,  i.e.  frcp=4640  Hz,  leads  to  the  Doppler  spectrum 
shown  in  Fig.  10. 
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Figure  9.  Doppler  spectrum  of  received  ground  clutter,  corresponding  to  two-way  array  pattern.  The  envelope 
represents  the  element  pattern. 

The  aliased  Doppler  spectrum  has  a  period  of  4640  Hz.  The  four  apparent  'periods’  observed 
within  the  4640  Hz  band  are  due  to  the  periodicity  of  the  grating  lobes,  and  we  note  that  their 
amplitudes  are  all  practically  equal. 


doppler  frequency  (kHz) 


Figure  10.  Aliased  Doppler  spectrum.  Pulse  repetition  frequency  f„p  chosen  such  that  the  4th  grating  lobe  in 
Fig.  8  aliases  onto  main  beam  (frep  =  4640  Hz). 


There  is  little  spreading  of  the  spectral  lines  associated  with  each  grating  lobe,  implying 
significant  periodicity  of  the  array  grating  lobe  Doppler  frequencies  over  the  beam-width  of  the 
satellite  array  pattern.  Also,  the  spectrum  between  the  grating  lobes  does  not  decrease 
monotonically  as  with  a  uniformly  illuminated  linear  array.  This  is  due  to  the  effective  amplitude 
modulation  imposed  by  the  two-dimensional  vertical  aperture. 

For  moving  target  indication  (MTI)  these  clutter  lines  are  simply  suppressed  by  an  N-pulse 
canceller  and  the  total  power  received  over  the  band  is  then  integrated.  The  canceller,  as  shown  in 
Fig.  5,  provides  an  output  signal 
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j2nfDnr 
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N 


Sou,='ZWnSn=SoYjWne 


n= I 
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where  wn  denote  weights,  and  the  successive  radar  pulses  sn  =  s0  exp(  j2  ji  fD  n  t ),  with  s0 
representing  the  amplitude  of  the  radar  pulse  and  2  n  fD  n  r  representing  the  phase  shift  due  to  the 
Doppler  frequency  fD  of  the  moving  target,  and  we  assume  x  =  l/frep. 

Choosing  the  usual  binomial  weights  for  the  canceller  leads  to 

/l  =  l 


or  a  normalized  canceller  power  transfer  function 


g(fD)  =  sin2(*  -l)(/c  fDr)  =  sin  2(W  x\nfD  /  frep)  .  (6) 

This  function  provides  nulls  at  fD=0  and  fD=freP,  i.e.  one  null  over  the  interval  (0,  fD).  A  larger 
number  of  evenly  spaced  nulls,  say  NT  nulls,  can  be  obtained  with  larger  time  delays  x  =  Nt  /frep. 
which  corresponds  to  using  every  NTth  pulse.  The  power  transfer  function  becomes 

g  (fD)  =  s\n2iN  -l)(7r  N  rfD  /  frep)  .  (7) 

Thus  we  have  two  parameters  to  suppress  the  clutter  spectral  lines:  the  number  of  weights,  N, 
which  controls  the  width  of  the  null  and  is  chosen  depending  on  the  aliased  grating  lobe 
bandwidth,  and  the  number  of  pulse  intervals,  NT ,  which  is  chosen  depending  on  the  number  of 
aliased  grating  lobes. 

We  demonstrate  the  effect  of  the  canceller  on  the  clutter  spectrum  in  the  example  above,  Fig.  10.  It 
has  four  lines  and  thus  we  set  Nt  =4.  Chosing  N=2  does  not  quite  suppress  the  clutter  lines  (see 
Fig.  1 1 ,  top),  but  with  N=3  we  reduce  the  clutter  residue  to  a  fairly  uniform  level  across  the  whole 
Doppler  band  (see  Fig.  1 1,  bottom).  Note  that  the  canceller  weights  were  chosen  deterministically 
using  only  knowledge  of  the  Doppler  frequencies  and  spectral  width  of  the  aliased  grating  lobes, 
but  without  knowledge  of  other  details  of  the  clutter  spectrum  at  the  input. 
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Figure  1 1.  Clutter  suppression  with  an  N-pulse  canceller;  “filter”  is  the  power  transfer  function  of  the 
canceller;  “clutter  ‘in’”  and  “clutter  ‘out’”  represent  clutter  power  spectrum  at  the  input  and  output  of  the 
canceller:  2-pulse  canceller  (N  =  2,  NT  =  4  )  at  top;  3-pulse  canceller  (N  =  3,  NT  =  4  )  at  bottom. 


4.  Comparison  with  Randomly  Thinned  Array 

To  determine  the  merit  of  our  approach,  we  compare  it  to  the  standard  approach,  represented  by  a 
19-element  random  sparse  array.  This  array  has  a  two-way  pattern  with  smeared  sidelobes  at  an 
average  level  l/(Neiemem)  =  -26  dB  and  operates  with  a  non-tailored  pulse  repetition  frequency. 

Taking  the  signal/clutter  ratio  (SCR)  as  the  performance  criterion,  we  compute  the  processing 
gain,  i.e.  the  ratio  of  SCRou,  and  SCRjn  at  the  output  and  input  of  the  N-pulse  cancellers  for 
different  numbers  N,  as  shown  in  Table  1 .  For  the  standard  approach  the  processing  gain  is 
practically  independent  of  the  number  of  processed  pulses,  since  the  sidelobe  pattern  and  thus  the 
clutter  Doppler  spectrum  are  essentially  flat.  Our  approach  gives  a  7  to  17  dB  advantage  since  the 
clutter  is  concentrated  around  a  few  narrow  Doppler  lines  and  can  be  selectively  suppressed. 


No.  of  pulses  N 

2 

3 

4 

Novel  approach 

16.3 

25.5 

26.7 

Standard  approach 

9.6 

9.7 

9.8 

A  dB 

6.7 

15.8 

16.9 

Table  1.  Processing  gain  SCRout/SCRjn  (dB) 
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5.  Performance  of  the  Rotating  Array 

The  above  examples  were  all  given  for  the  elliptical  array  at  an  array  rotation  angle  (pa=  0°,  i.e. 
the  situation  shown  in  Fig.  6b  with  five  elements  on  the  major  axes  of  the  ellipses.  This 
represents  the  best  periodic  case.  However,  during  one  Earth  orbit,  the  array  rotates  360°  around 
its  center  and  changes  its  triangular  grid  and  thus  its  grating  lobe  structure.  Therefore,  we  also 
evaluated  performance  for  other  array  rotation  angles  tpa,  over  the  range  0  <  tpa  <  30°  in  2.5° 
increments,  after  which  the  array  configuration  repeats  due  to  symmetry. 

We  found  that  we  can  always  select  a  basic  pulse  repetition  frequency  that  leads  to  signal-to- 
clutter  (S/C)  improvement  factors  around  25  dB.  Fig.  12  shows  these  improvement  factors  and  the 
corresponding  PRFs  that  maximize  signal  to  clutter.  Here  S/C  improvement  was  defined  as  S/C  at 
the  filter  output  relative  to  S/C  at  the  single  receive  element,  assuming  target  Doppler  frequency 
to  be  uniformly  distributed  between  0  Hz  and  PRF,  and  the  filter  being  a  binomial,  5-pulse 
canceller.  Multiples  of  these  basic  PRFs  give  the  same  result  and  thus  our  approach  is  robust  with 
respect  to  both  array  rotation  angle  and  PRF. 


Figure  12.  Signal/Clutter  improvement  factor  vs.  array  rotation  angle  (top),  and  corresponding  pulse 
repetition  frequency  (bottom). 
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6.  Optimization  of  Separable  Space-Time  Processing 

Our  proposed  approach  for  separable  space-time  processing  has  M  receive  element  weights  am 
and  N  filter  weights  wn,  (see  Fig.  5a),  which  in  the  examples  above  were  chosen  corresponding  to 
a  maximum  gain  array,  (i.e.  am=l),  and  to  a  classical  N-pulse  canceller,  (i.e.  wn=binomial 
coefficients),  respectively.  However,  these  weights  can  be  chosen  so  as  to  maximize  the  S/C 
improvement  factor.  Optimizing  either  the  weights  am  or  the  weights  vv„  leads  to  a  ratio  of 
quadratic  forms  in  terms  of  these  weights,  which  has  a  well-known  maximum  solution  [22]. 


We  have  performed  such  computations  for  various  filter  orders  and  the  array  rotation  angles, 
<pa=0°  and  cpa=5°,  which  represent  best  and  worst  cases,  respectively.  The  results  are  shown  in  Fig. 
13.  We  note  that  array  optimization  is  considerably  more  effective  than  filter  optimization.  Also, 
apparently  the  N-pulse  canceller  matches  optimum  weighting  for  low  filter  orders  (2-4  pulses)  but 
degrades  for  higher  orders.  This  is  because  the  deterministic  weights  place  nulls  only  over  the 
dominant  Doppler  lines.  Once  these  first  order  contributors  have  been  suppressed  to  a  uniform 
level,  additional  nulls  at  these  locations  will  not  improve  the  processing  gain.  Optimized  weights, 
however,  continue  to  null  second-order  contributors,  so  performance  improves,  although  slowly. 


S/C  Improvement  Factor  vs  Filter  Order 


Figure  13.SignaI/Clutter  improvement  factor  vs.  filter  order  for  various  alternatives  of  optimization. 


A  joint  optimization  of  the  array  and  filter  weights  can  be  achieved  in  an  iterative  fashion  by 
alternating  optimizing  the  am  and  wn.  Convergence  is  rapid  and,  although  not  guaranteeing  a 
global  maximum,  leads  to  a  maximum  improvement  factor  attainable  with  separable  space-time 
processing.  The  results  shown  in  Fig.  13  represent  only  slight  improvement  over  the  optimized 
array  with  a  simple  N-pulse  canceller. 


7.  Comparison  with  General  Space-Time  Processing 

Finally,  for  reference,  we  compare  separable  processing  with  the  ultimate  case  of  space-time 
optimal  processing  (STOP),  where  each  array  element  is  followed  by  an  individual  Doppler  filter. 
In  contrast  to  separable  processing  which  employs  M+N  degrees  of  freedom,  STOP  employs 
MxN  degrees  of  freedom.  Performing  an  analysis  as  before,  but  now  based  on  the  full  space-time 
clutter  covariance  matrix,  leads  to  the  results  shown  also  in  Fig.  13.  The  additional  degrees  of 
freedom  clearly  give  superior  improvement  factors.  When  the  number  of  pulses  exceeds  3,  the 
curves  are  noise-limited. 
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Our  separable  approach  is  critically  based  on  a  periodic  array  and  a  ‘tuned’  PRF.  An  interesting 
question  therefore  is  whether  the  flexibility  of  general  space-time  processing  can  relieve  this 
requirement  for  array  periodicity. 

To  explore  this  question,  we  added  random  errors  to  the  element  location,  such  that  they  were 
distributed  uniformly  within  a  square  5  m  box  in  the  orbital  plane  and  centered  at  the  desired 
points.  Assuming  that  the  errors  were  known  by  independent  measurements  we  then  recomputed 
the  S/C  improvement  factors  for  these  arrays  with  perturbed  periodicity.  The  results  for  two 
sample  cases  are  shown  in  Fig.  14.  Clearly  even  the  general  STOP  solution  benefits  from  a 
periodic  array  lattice. 


Figure  14.  Comparison  of  Signal/Clutter  improvement  factors  achievable  with  a  periodic  array  and  array  with 
periodicity  perturbed  by  random  but  known  errors  (two  realizations). 


8.  Performance  of  the  Planar  Array  under  Non-Ideal  Conditions 

In  the  previous  sections  we  have  analyzed  a  vertical  planar  periodic  array  at  one  single,  carefully 
chosen  pulse  repetition  frequency,  and  have  demonstrated  this  to  be  an  ideal  case  with  optimum 
clutter  suppression.  However,  controlling  the  satellites  to  maintain  the  exact  positions  of  a 
periodic  grid  may  not  be  possible.  Therefore  we  now  analyze  how  performance  degrades,  first 
with  array  pitch  and  yaw  and  second,  with  random  but  known  element  position  errors.  We  do  this 
as  a  function  of  pulse  repetition  frequency,  since  this  is  our  main  free  variable. 

8.1.  The  Effects  of  Array  Pitch  and_  Yaw 

We  consider  the  same  19-satellite  cluster  as  before,  which  forms  a  planar  vertical  array  as  shown 
in  Fig.  6.  For  this  array,  ‘pitch’  implies  that  the  velocity  vector  differs  from  the  array  axis 
direction  but  lies  in  the  plane  of  the  array.  It  is  caused  by  the  array  rotating  around  its  center 
during  orbit.  ‘Yaw’  implies  a  velocity  component  orthogonal  to  the  plane  of  the  array.  It  is  caused 
by  the  Earth’s  rotation. 
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Figure  15.Improvement  factor  vs.  pulse  repetition  frequency  for  the  19-element  array  with  0  deg.  pitch,  0  deg 
yaw  (top),  and  5  deg.  pitch,  0  deg.  yaw  (bottom). 

For  reference  we  show  in  Fig.  15  (top)  the  improvement  factor  (IF)  vs.  pulse  repetition  frequency 
for  0  deg.  pitch  and  0  deg.  yaw.  This  is  the  optimal  case.  We  note  that  there  are  many  alternative 
prf  values  which  give  quite  large  IF  values.  However,  these  values  sometimes  have  an  extremely 
narrow  bandwidth,  e.g.  for  the  peak  at  prf=198  Hz,  a  0.005  Hz  change  makes  a  10  dB  difference. 
The  reason  is  that  the  peaks  essentially  represent  the  difference  between  two  large  numbers,  the 
clutter  returns  from  two  successive  pulses.  Since  we  cannot  sample  the  curve  densely  enough 
these  IF  maxima  are  not  always  properly  resolved,  and  in  fact,  the  true  IF  maximum  =73  dB  at 
198.027  Hz. 

The  lowest  prf  for  which  an  IF  maximum  occurs  corresponds  to  the  situation  when  just  the  two 
edge  elements  of  the  array  form  a  DPCA  pair,  i.e.  at  prf=2  VR/Ad  =  148.5  Hz,  for  a  satellite 
velocity  V«  =  7425  m/s  and  element  distance  Ad  =  aperture  size  =  100  m,  which  yields  an  IF  =73 
dB.  At  higher  prfs  some  still  larger  IF  values  occur  due  to  the  simultaneous  formation  of  several 
DPCA  pairs. 

Fig.  15(  bottom)  shows  IF  vs.  prf  when  the  array  has  a  pitch  angle  of  5  deg.,  which  represents  one 
of  the  more  difficult  cases.  We  note  that  now  there  are  many  more  IF  maxima  and  also  that  the 
minima  are  less  deep.  This  is  due  to  the  fact  that  the  rotated  array  presents  a  larger  set  of 
interelement  spacings  and  thus  a  larger  set  of  DPCA  pairs. 
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Figure  16.  Improvement  factor  vs.  pulse  repetition  frequency  for  the  19-element  array  with  0  deg.  pitch,  and 
yaw  corresponding  to  Earth  rotation. 


Fig.  16  shows  the  effect  of  yaw  due  to  Earth  rotation.  Apparently,  this  has  very  little  effect  on  the 
IF  curve.  Presumably  the  reason  is  that  the  array  elements  are  highly  directive  and  thus  the  returns 
come  essentially  from  the  same  direction,  all  with  the  same,  relatively  small,  extra  Doppler  shift 
caused  by  the  Earth  rotation.  The  last  observation  that  Earth  rotation  has  only  a  minor  effect  holds 
also  for  a  5  deg.  pitch  angle. 

8.2.  Effects,  of  Random  Position  Errors^ 

In  this  section  we  investigate  array  performance  degradation  when  the  array  (satellite  cluster)  can 
not  be  controlled  to  maintain  a  perfectly  periodic  grid.  To  model  this  situation  we  distribute  the 
array  elements  randomly  within  square  boxes,  which  are  centered  at  the  periodic  grid  points  and 
contained  in  the  single  orbital  plane  of  the  array.  Assuming  the  position  errors  are  known,  we 
again  compute  the  improvement  factor  as  a  function  of  pulse  repetition  frequency  for  various 
error  box 
sizes. 

Fig.  17  shows  the  improvement  factor  vs.  pulse  repetition  frequency  for  our  19-element  array 
with  0  pitch  and  0  yaw  for  error  boxes  1 , 3,  and  5  m  on  the  side,  for  one  random  but  known 
realization.  The  case  of  0  m  box  size,  no  position  errors,  shown  in  Fig.  15  (top)  above  serves  as  a 
reference.  With  increasing  position  errors,  the  number  of  IF  maxima  increases  while  their  peak 
values  decrease.  This  is  again  because  there  is  a  richer  set  of  element  spacings,  but,  at  any  one 
frequency,  only  a  small  number  of  elements  contribute  to  form  DPCA  pairs. 

Finally,  Fig.  18  shows  IF  results  which  include  the  effects  of  a  5  deg.  pitch  angle  and  yaw  due  to 
Earth  rotation.  Clearly  the  increased  randomness,  caused  by  the  array  rotation  and  the  increasing 
error  boxes,  severely  limits  the  attainable  improvement  factor.  It  is  difficult  to  give  an  accurate 
estimate  for  this  degradation  due  to  the  large  dynamic  range  of  the  IF  curves.  However,  it  appears 
that  random  errors  associated  with  the  5  m  box  cause  an  IF  loss  on  the  order  of  10  dB  compared 
with  the  ideal  case  with  no  position  errors. 
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9.  Conclusions 

We  have  proposed  a  novel  two-dimensional  periodic  array  geometry  for  TechSat21  and  an 
associated  technique  for  pattern  synthesis  in  angle-frequency  space.  It  exploits  the  double 
periodicities  of  the  grating  lobes  in  the  angular  domain  and  of  the  radar  pulses  in  the  frequency 
domain,  and  allows  substantial  gains  in  clutter  suppression. 

We  have  considered  a  ‘separable’  system  architecture  for  TechSat21,  where  a  planar  M-element 
array  is  followed  by  a  single  N-pulse  Doppler  filter,  and  have  evaluated  some  processing 
alternatives  for  their  effectiveness  in  clutter  suppression.  An  attractive  feature  of  this  architecture 
is  that  it  allows  simple,  multi-pulse  MTI  processing  and  adaptive  antenna  pattern  control.  The 
limited  number  of  degrees  of  freedom,  M+N,  require  a  relatively  low  computational  load. 

General  space-time  processing  offers  substantially  improved  signal/clutter  ratios.  However,  it 
employs  a  much  larger  number  of  degrees  of  freedom,  MxN,  and  may  be  associated  with  a 
computational  load  that  is  difficult  to  realize.  Synthesis  techniques  for  thinned  apertures  on 
transmit  and  receive,  and  the  trade  between  performance  and  processing  complexity  merit  further 
study. 

Finally,  taking  the  ideal  periodic  array  as  a  reference  case,  we  have  analyzed  the  performance 
degradation  with  array  pitch  (satellite  cluster  rotation)  and  yaw  (Earth  rotation)  and  known 
random  element  position  errors.  We  found  that  the  signal/interference  improvement  factor  is  an 
extremely  sensitive  function  of  the  pulse  repetition  frequency  and  that  it  varies  over  a  large 
dynamic  range.  Achievable  improvement  factor  maxima  are  moderately  degraded  by  array  pitch 
and  strongly  degraded  with  increasing  position  errors.  Yaw  effects  appear  to  be  negligible.  In 
general,  with  increased  randomness,  fewer  elements  contribute  to  form  DPCA  pairs  so  that 
satellite  usage  efficiency  (and  SNR)  is  degraded. 

Thus,  a  periodic  array  lattice  allows  for  significantly  better  clutter  suppression  than  a  non-periodic 
lattice,  independent  of  the  particular  space-time  processing  scheme. 
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Figure  17.  The  improvement  factor  vs.  pulse  repetition  frequency  for  random  position  errors  within  boxes  1,  3, 
and  5  m  on  the  side  (0  deg.  pitch,  0  deg.  yaw). 
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Figure  18.  Improvement  factor  vs.  pulse  repetition  frequency  for  the  case  of  array  pitch  =5°,  yaw  according  to 
Rarth  rotation,  and  increasing  element  position  errors. 
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Appendix  A 

Comparison  of  Target  Detection  for  Two  Different  Radar  Modes 

We  show  that  the  target  detection  performance  against  receiver  noise  of  an  array  radar  where  each 
element  transmits  its  own  orthogonal  signal  can  be  equivalent  to  that  of  a  conventional  radar  with 
a  single  coherent  transmit  signal.  To  this  end  we  compare  the  total  time  required  by  the  two 
systems  to  survey  a  given  angular  sector  Qs ,  at  a  given  signal-to-noise  ratio. 

Both  systems  are  assumed  to  have  an  array  of  N  elements  with  per  element  transmit  power  PQ, 
transmit  gain  G,,  receive  gain  Gr,  and  thermal  receiver  noise  no.  The  element  pattern  beamwidth  is 
Q0  and  the  array  beamwidth  is  £2a ,  with  Qa «  Q0  (see  Fig.  A 1 ). 

We  consider  a  surveillance  sector  equal  to  the  element  beamwidth,  Qs  =Q0,  which  is  motivated  as 
follows.  The  satellite  cluster  constitutes  a  phased  array  where  each  element  is  one  of  the  phased 
arrays  on  the  individual  satellites.  Thus,  in  contrast  to  a  conventional  array  whose  element 
patterns  are  fixed  in  space,  in  the  cluster  array  the  element  patterns  can  be  redirected.  However, 
once  these  directions  have  been  set  and  the  elements  illuminate  a  fixed  sector  of  width  Q«,  then  the 
cluster  array  is  equivalent  to  a  conventional  array.  Surveillance  is  then  simply  performed  by 
scanning  with  the  narrow  cluster  array  beam  of  width  £2a  over  the  element  pattern  sector  Q0. 
Surveillance  over  a  sector  larger  than  the  element  (satellite)  pattern  sector  £20  is  accomplished  by 
step  scanning  the  element  pattern  to  contiguous  sectors  of  width 

1)  Conventional  Radar  Mode.  Here  the  array  elements  (satellites)  transmit  coherent  pulsed  signals 
all  on  the  same  carrier  frequency.  At  the  target,  the  illumination  intensity  is  (NP0)NGt,  and  for  a 
fixed  radar  cross  section,  the  power  received  is  oc  N'  P„G,  Gr. 

After  coherent  integration  the  receiver  noise  is  n„  relative  to  the  signal  power  received.  Thus  the 
signal-to-noise  ratio  SNR  °c  N  P0GtGr/  n0,  and,  assuming  equal  power  and  dwell  time  in  each 
beam  position,  the  total  surveillance  time  is  «=  £2,,/  Qa. 

2)  Multi-Signal  Mode.  In  this  case  the  elements  (satellites)  simultaneously  transmit  orthogonal 
signals,  e.g.  pulses  on  different  carrier  frequencies  fi  to  f^  It  is  assumed  that  the  relative  locations 
and  pulse  timings  are  so  precisely  known  for  each  element  that,  after  coherent  demodulation,  the 
complex  video  pulses  corresponding  to  each  orthogonal  signal  received  at  one  element  from  one 
target  direction  can  be  combined  coherently.  Furthermore,  it  is  assumed  that  the  spread  of  carrier 
frequencies  fi  to  fN  is  not  so  large  as  to  decorrelate  the  complex  signal  scattered  by  the  target.  In 
effect,  this  is  coherent  operation  on  transmit,  and  generates  an  effective  illumination  function 
equal  to  the  array  pattern.  The  following  analysis  assumes  pulses  transmitted  on  different  carrier 
frequencies,  fi  to  fN,  but  applies  as  well  for  general  orthogonal  signals  such  as  might  be  generated 
by  orthogonal  phase  coding  of  the  transmitted  pulses. 

Element  1  will  illuminate  the  sector  Q„  with  intensity  «  P0G,  at  carrier  frequency  fi.  The 
received  signals  at  each  of  the  N  array  elements  from  the  element  1  transmission  can  be  combined 
coherently  to  form  a  beam  of  width  Qa  with  gain  NGr  and  a  SNR  «  NP„G|Gr/  no  .  However,  since 
the  entire  sector  £2C  is  illuminated,  we  can  form  a  set  of  N  contiguous  receive  beams  within  this 
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sector  via  digital  beam  forming.  This  process  can  be  repeated  at  the  frequencies  f|  to  fN,  with  one 
frequency  transmitted  from  each  of  the  N  elements,  thus  allowing  a  total  of  N2  independent 
receive  beams.  For  the  sparse  array,  the  number  of  independent  beam  positions  within  the  sector 
£20  is  not  exactly  N  (~  20)  but  rather  on  the  order  of  50,  based  on  the  size  of  the  sparse  array  and 
the  along-track  aperture  of  the  satellite  array.  This  difference  will  be  neglected  here. 

These  beams  can  be  employed  in  different  ways: 

a)  Form  a  cluster  of  N2  contiguous  beams.  For  the  desired  surveillance  sector  this  leads  to  a 
surveillance  time  °e  Q0/  N"  £2a  at  a  SNR  «=  NP0G,Gr/  n0.  In  comparison  with  the  conventional 
radar  mode,  both  the  surveillance  time  and  the  signal  to  noise  ratio  are  reduced  by  N2. 

b)  Form,  at  each  of  the  N  frequencies,  the  same  overlapping  cluster  of  N  contiguous  receive 
beams,  and  add  the  N  frequency  signals  at  each  beam  output  coherently.  When  compared  to 
option  (a),  this  leads  to  an  increased  surveillance  time  °cQ0/NOa  at  a  proportionally  increased 
SNR  «  N‘P0G,Gr/  nc.  Thus,  in  comparison  with  the  conventional  radar  mode,  the  surveillance 
time  and  the  signal-to-noise  ratio  are  reduced  by  N. 

c)  Finally,  integrating  coherently  N  time  samples  from  the  above  N  beam  clusters  gives  a 
surveillance  time  Q0/  Qa  at  a  SNR  «  N3  P„G,Gr/  n0  . 

This  last  alternative  shows  that  the  conventional  radar  mode  and  the  multi-signal  mode  give  the 
same  basic  detection  performance  against  receiver  noise.  However,  the  latter  mode  allows  a  trade 
of  surveillance  rate  vs.  SNR  and  longer  integration  times  for  Doppler  processing  to  suppress 
clutter. 


Q„ 


-► 

angle 


Figure  Al.  The  element  beamwidth  is  much  wider  than  the  array  beamwidth  Qa.  Surveillance  of  a  radar 
sector  a  =  no  requires  £VQa  beams. 
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Appendix  B 

The  Doppler  Spectrum  for  Moving  Targets  and  Clutter 


B 1 .  Aperture  and  Integration  Time  Requirements  for  Ground  Moving  Target  Detection 

Here  we  derive  requirements  on  antenna  aperture  and  coherent  integration  time  to  detect  slowly 
moving  ground  targets. 

We  assume  the  radar  has  a  horizontal  velocity  vradar  and  illuminates  one  circular  range  cell  at  a 
depression  angle  8,  as  shown  in  Fig.  B 1 .  The  vector  vtargc,  denotes  the  target  velocity  relative  to 
ground,  and  v  is  the  angle  between  the  target  velocity  and  the  direction  R  to  the  radar.  From  the 
same  azimuth  angle  (p  we  thus  get  two  different  Doppler  returns  due  to  the  target  motion,  as 
illustrated  in  the  Angle-Doppler  diagram  in  Fig.  B2. 


Figure  Bl.  Illustration  of  radar  and  moving  target  geometry. 


The  Doppler  frequency  from  ground  clutter  at  azimuth  angle  cp  is 


f  ^  hraJar  \  R 
J  clutter  n 


leading  to 


where 


f clutter  fc max  V 

fCmx  =_^cos  8. 

J  t  max 


(Bl) 

(B2) 
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The  Doppler  frequency  for  the  target  is  similarly 


•f 'larger 


leading  to 


^  radar  +V,arRe<)-R 


2v„ 


f  target  =  fcluner  +^y^COS/9  . 


(B4) 


(B5) 


Defining  the  angle  a  between  the  vector  v,argci  and  the  tangent  to  the  clutter  cell  we  find  that 

cos  (i  =  cos  <?si  n  a  (B6) 

and  thus  the  Doppler  differential  A/iarge,  between  the  clutter  and  the  target  is 


^f target  ft arget  ft flutter 


 ~ ''  target 


cos^sinor. 


(B7) 
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Figure  B2.  Target  return  and  ground  clutter  return  in  Angle-Doppler  space. 
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We  can  now  estimate  the  minimum  aperture  size  and  integration  time  as  illustrated  in  Fig.  B3. 


Figure  B3.  Illustrating  the  derivation  of  minimum  aperture  size  and  coherent  integration  time. 


For  a  given  Doppler  differential  A  /target  determined  by  the  minimum  detectable  target  speed  we 
first  require  that  the  pass-band  clutter  be  in  the  first  sidelobe  of  the  antenna  pattern,  i.e. 


A  (p- 


3  A 

2  \mJ 


(B8) 


where  Aproj  =  A  costp  is  the  projected  antenna  aperture  and  A  is  the  aperture  along  the  direction 
of  vradar.  From  (B2)  above  we  obtain 


/o*CO.S«> 

A  (p 


(B9) 


and  setting 


clutter  =A/„ 


larger 


(BIO) 


we  obtain  from  (B3)  and  (B7) 


A  = 


3  VruJar 

2  »W,sin  a 


A. 


(Bll) 
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Similarly  we  require  that  the  main-beam  clutter  be  in  the  first  sideband  of  the  Doppler  filter,  i.e. 


44 


target 


3J_ 

2  T 


(B12) 


where  T  is  the  coherent  integration  time. 
Using  (B7)  and  (B8)  leads  to 


T 


radar 


4  v,arg',  cos  S sin  a 


A. 


(B 13) 


In  (B 1 1 )  and  (B 1 3)  the  dependence  on  the  angles  a  and  5  represent  the  unavoidable  angular 
dependence  of  any  Doppler  detection  scheme.  Ignoring  these  factors  we  thus  find  for  the 
minimum  antenna  aperture  size 


2  V, 


target 


and  for  the  coherent  integration  time 


<2 

l )>£  2™^ .  X. 


4  v, 


target 


(B14) 


(B15) 


B2.  The  Doppler  Spectrum  of  Ground  Clutter  Seen  from  a  Satellite  in  Earth  Orbit 

The  Doppler  frequency  of  the  ground  reflection  as  seen  from  a  satellite  in  a  circular  Earth  orbit  is 

derived  in  [18]  and  repeated  here  for  reference. 

With  the  following  notation 

i  =  angle  between  orbital  and  equatorial  planes,  see  Fig.  B4 

v  =  angular  position  of  the  satellite  in  the  orbital  plane,  measured  from  the  equator,  see  Fig.  B4 
Rc  =  radius  of  the  Earth  =  6378  km 

hs=  height  of  the  satellite  above  the  surface  of  the  Earth  ~  850  km  for  our  examples 

Pe  =  Earth  gravitational  constant  =  3.986xlOIJ  m3  /sec2 

(0e  =  angular  rate  of  Earth  =  7.29x  10'5  rad/sec 

(p  =  azimuth  angle  with  respect  to  array  broadside,  see  Fig.  B4 

5  =  depression  angle  to  the  ground  reflection  point,  see  Fig.  B4 

dR/dt  =  line  of  sight  distance  change  with  time 

X  =  wavelength  of  radar  signal 

the  Doppler  frequency  of  the  clutter  is 
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(B16) 


=  _2dR_2 
D~  A  dt~  A 


R+h 


sin  q>  cos  8 -  0)t  ( Rr  +  hs  )(sin  <p  cos  i  -  cos  ^  sin  i  cos  V)  cos  S 


Figure  B4.  Global  coordinate  system  and  satellite-based  coordinate  system  used  in  the  computation  of  the 
clutter  Doppler  frequency  (rs  and  vs  are  unit  vectors  along  the  radial  direction  to  the  satellite  and  along  the 
satellite  velocity,  respectively,  and  n5=r#\v,). 
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Appendix  C 

Optimization  of  the  Array  Element  Weights 

In  this  appendix  we  derive  the  optimized  array  weights  {am},  such  that 


P„„,  (a>  £'  =  signal) 

Pout  (a’ E  =  clutter)  +  Pul  (noise) 


max.. 


(Cl) 


where  Pout  is  the  Doppler  filter  output  power,  E1  denotes  the  incident  field  on  the  array,  and  the 
Doppler  filter  weights  { wn}are  assumed  given  (see  Fig.  Cl). 


Ein 


Figure  Cl.  Array  beam  former  with  weights  am  followed  by  Doppler  filter  with  weights  wn 

We  consider  a  plane  wave  E‘(<p)ejW'^k  r)  incident  on  an  array  with  elements  m  located  at  rm . 
Here  f=fo+fo  denotes  the  sum  of  the  basic  carrier  frequency  fo  and  the  Doppler  frequency  fo  . 
The  direction  of  the  wave  vector  k  denotes  the  direction  9  of  the  wave  *,  and  its  magnitude 
|  k  |=  2 nflc,  where  c  is  the  speed  of  light.  After  demodulation  (which  essentially  strips  off  the 

carrier  frequency)  the  received  signal  at  element  m  is  E‘((p)  eimfn'~v  Tm)e0(<p)  ,  where  e«(cp) 

denotes  the  element  pattern  with  phase  center  at  the  origin  (assumed  equal  for  all  elements).  The 
array  output  voltage  thus  becomes 

v0(/,/D)  =  E'(<P)  ej<0n‘  e0(<p)  ,  (C2) 

m 

where  am  are  the  element  weights. 


*)  The  use  of  a  single  angular  variable  is  justified  since  a  range  bin  maps  onto  a  very  narrow  sector  in  the  array 
pattern,  as  discussed  in  the  context  of  Fig.  8.  Thus  the  total  received  clutter  power  can  be  obtained  by  integration  over 
a  single  angular  variable  (the  azimuth  angle  cp)  along  the  sector,  while  assuming  the  power  density  being  constant  in 
the  orthogonal  direction. 
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In  all  following  computations,  we  assume  k(/)  (rm  -r.)  -  k(/0)  (rm-r„),  since  the  Doppler 
frequencies  and  range  differences  are  sufficiently  small. 


Defining  column  vectors  a  and  <1>  with  components  am  and  ,  respectively,  we  can  write  the 
amplitude 


va(fD)  =  E‘e0(<p)a'  <I>(k)  .  (C3) 

The  Doppler  filter  has  the  transfer  function  w(fn )  so  that  the  output  voltage 

Ku,  =  w(fD)  va(v.fD)  =  E‘  w(fD)  e0(<p)  a'  <D(k) .  (C4) 


For  the  target  signal  we  assume  the  incident  field  ETe~Jkr'r  to  come  from  a  single  target  direction 

<pT  with  amplitude  E'  uniformly  distributed  over  the  frequency  interval  (0,  frcp).  This  leads  to  a 
total  output  power 


^=a'<2a\ 


where  the  elements  qmn  of  Q  are 


=  I ET 


k<M)  I2  }  K/0) 


I2  K-'J 


dfD 


Recalling  that  kr  •  (rm -r„) |»kr(/0)-(rM -rn) |  weobtain 


qmn  ~  constant xe  y1,r(/,’>  <r"-r") . 

Thus  the  total  target  signal  power 

PT  =a'  Qa  oc  a'<I»(k7.(/0))<^(k7.(/0))ta*. 


(C5) 


(C6a) 


(C6b) 


(C7) 


(C8) 
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The  incident  field  due  to  clutter  is  spread  over  both  angle  (p  and  Doppler  frequency//? .  However, 
we  note  that  there  is  a  one-to-one  correspondence  between  (p  and  fD,  <p=(p(fD ),  and  furthermore, 
the  clutter  signal  is  incoherent  over  <p. 

The  power  contribution  A P(  received  from  a  sector  A<p  is  proportional  to  the  sector  size,  the 
clutter  cross  section  per  angle,  <70,  (assumed  constant),  the  relative  illumination  intensity  given  by 

the  transmit  power  pattern  |  Ft  \ 2  and  the  receive  power  pattern  |  Fr  |2 ,  i.e. 


APC  -e^F,^)?  A#\Fr{(p)\\ 


Thus  the  received  clutter  power  spectrum 


Sc 


APC 

4/o 


oilW0)PIW®)P-^ 

dfD 


or,  with 


<*c  (  /d  ) — 


^oI^(/d)|2 


(C9) 


(CIO) 


(Cl  1) 


SC(fo) 


^(/o)l^/0)l2=^(/fl)ko(/D)ZvjVr'f.  (Cl  2) 


Using  similar  vectors  a  and  O  as  before  we  obtain 

Sc(fD)  ~  <rc(fo)\e0(fo)\2  *  <WfDY a  .  (C13) 


The  matrix  A/0(/o)  =  <&(fD  )<!>(/„ )+  is  computed  numerically,  as  a  function  of  frequency,  by 

computing  for  each  element  location  rm  a  set  of  values  e~jk<fDm)  r"  ,  where  the  sample  frequencies 
fDn  extend  over  the  entire  Doppler  spectrum.  (  The  wave  number  k(fDn)  has  magnitude  2n/Xo 
and  direction  (p  =  (p(fon  ))• 

The  aliased  Doppler  spectrum  Suc/(tt(  is  obtained  as 


SL  =  £s  c</„ +///„)= 


(Cl  4) 


=  a 


£(TC(/d  +/lfrep)\e0(fD  +MfreP)\2  MMd  +Vfr<p) 


which  leads  to  the  total  clutter  output  power 
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(Cl  5) 


Jrep  oo 

rf„=  JN/0)I!  s;lu,(f„)dfD  =  l\w(fD) f  sc(fD)df0, 


since  w(fD  +///  )  =  w(fD)  for  all  //.  Thus 


=a'Ca* 


where  the  clutter  matrix 


oo 

C=  f  I  w(/)|2  <rc(/D)|e0(/i>)|2  dfD  . 


Finally,  adding  thermal  noise  PN  =a'  nj  a*  to  the  clutter  gives  the  signal 

*£  +  /£  a'(C  +  n0/)a 

This  ratio  of  quadratic  forms  is  maximized  by  the  weight  vector 
a  =  a[(C  +  «„/)-'  O(/0)]*. 


(Cl  6a) 

(Cl  6b) 

to-clutter+noise  ratio 

(Cl  7) 

(Cl  8) 
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Appendix  D 

Doppler  Filter  Optimization 

For  an  array  followed  by  a  Doppler  filter,  as  shown  in  Fig.  Dl,  we  determine  the  filter  weights  wn 
such  that  the  signal-to-clutter  ratio  is  maximized  at  the  output  of  the  Doppler  filter.  We  use  the 
following  notation 

vm(/o’0  =  complex  voltage  at  filter  input 
v»w(/d»0  =  voltage  at  filter  output 
wn  =  filter  weights 

f 0  =  Doppler  frequency 

t  =  time  between  pulse  samples,  normally  x  =l/frep,  i.e.  every  consecutive  pulse 

is  used,  however,  setting  x  =Nt/  frep  allows  summation  over  every  NTth 
pulse  only. 


Figure  D  1.  Array  with  Doppler  filter. 

Sampling  the  demodulated  input  signal  at  times  t  =  (n-1)  x  ,  n=l  ,2,...,  N,  and  delaying  and 
summing  these  pulses  we  obtain 

vou, (/o )  =  Mfo )  Z  vv«  ej2*in~')foT  =  w'  v,.„  (D 1  a) 

n=l 

where 

v'n  =  A(/D)(  le**" . 

w 1  =(wt,w2,...,wN).  (Dlb) 

The  output  power  spectrum  Sout  is  proportional  to  the  square  of  voul,  i.e. 

SoU,  =\MfD)\2  w' M  w'  (D2a) 

where  the  matrix 
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M(fD)  = 


(D2b) 


r\  e-i**fo*'  g-j2*(N-DfnrSl 

eilKfnT  l 


J2x(N-\)fnT  , 

ye  .  1  J 


Now,  assuming  that  the  target  Doppler  frequency  is  uniformly  distributed  over  the  interval  (0, 
hep),  he.  A  (fD)  =  constant,  and,  since  r  =  Nr  /  frrp ,  we  obtain  for  the  total  target  signal  power 


PT  = 


'  rrp 

|  ST„u,(fD)dfD  =\A'\  w' 


w'w\  (D3) 


Similarly,  we  obtain  for  the  total  clutter  power 

f rrp  fret 


P»,  =  J  S^ifoWo  =  w'  J  |ac(/d)|2  M(/0)J/d  w* 


(D4) 


where  | A(  (fD) j  denotes  the  aliased  clutter  power  spectrum  at  the  filter  input,  as  given  by  (C14). 

frrp  ^ 

The  coefficients  in  the  integral  matrix  A/  .  =  J  Ac(/D)|"  M(fD)dfD  are  of  the  type 


Jrep 

cm  =  J  |ac(/0)| 


(D5) 


and  have  to  be  evaluated  numerically  (they  correspond  to  the  Fourier-series  coefficients  of 
| A*  (/„) |  ).  With  this  notation  we  get 


A*c  =  /, 


rep 


C0  C-\  C-2 
C\  CQ  C_, 

c2  c,  c0 


(D6) 


which  is  Hermetian,  since  c_m  =  cm  .  The  total  output  clutter  power 


pL  =  w'  Mc  w‘ 


(D7) 


We  want  to  maximize  the  signal-to-clutter  ratio,  i.e.  we  seek  a  weight  vector  w  such  that 
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pL  »v 

=—  (D8a> 


or,  equivalently, 


1 

— - —  =  maximum 

w'  Mc  w 


(D8b) 


where  w  denotes  a  weight  vector  of  unit  magnitude. 

The  solution  to  this  maximization  problem  is  given  in  terms  of  the  eigenvectors  e(1  and 
eigenvalues  Xn  of  Mr ,  [22],  and  is 


w 


PT 

out 

pc  , 
rom  y 


(D9) 


where  ‘min’  denotes  the  minimum  eigenvalue  and  the  corresponding  eigenvector. 
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Appendix  E 

Optimization  of  Space-Time  Filter 


In  this  appendix  we  derive  the  optimized  weights  { wmn}  of  a  space-time  Filter,  such  that 


p„u, (K„„K£'  =  target  signal) 
p,m, ({wmn}’E>  =  clutter )  +  P  (noise) 


=  max. 


(El) 


where  Pout  is  the  filter  output  power  and  E'  denotes  the  incident  field  on  the  array  (see  Fig.  El). 


^ - 1 

t\ 

►fwiTI— 


r 


^ — l 

t\ 


T 


Wmn 


out 


Figure  E  1.  Array  with  full  space-time  Filter. 


We  consider  a  plane  wave  E'e  Jk  r  incident  on  an  array  with  element  m  located  at  rm ,  as  before. 
After  demodulation  the  received  voltage  is  at  element  m,  time  delay  n 

v™  =v„  (/-(*-!  )T)  =  Eiem«p)eJUM,-(n-')T)  (E2) 


where  the  element  pattern  em(<p)  =  e0(tp)  e  Jk  r”  and  e0(<p)  denotes  a  common  element  factor. 
The  voltage  output  from  the  space-time  filter 

v»u,  =Z«,™V™  =  W'V  (E3) 

m,n 

where  the  MxN-dimensional  vectors  w  and  v  have  the  components  w„  =  w(m_l)A,+n  and 
=  v(m.1)W+n ,  respectively. 
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The  total  power  received 


frrp 

=  w'  |  v  vV/„  w' 

0 


(E4a) 


or 

^,u,  =w'M  w\  (E4b) 

where  the  elements  of  the  covariance  matrix  M  are 

W4fo  ■  (E4c) 

0 


For  the  target  we  assume  the  signal  to  come  from  angle  <pr  and  be  uniformly  distributed  over  the 
Doppler  range  (0,  frep),  as  before.  This  leads  to  a  target  covariance  matrix  M r  with  elements 


0 

*  |£reo(0V)|2  e~ikT''hHtmW~tm{v))  J  e-i2’'<nw-n<v»fTdfD  = 


(E5) 


—  ^o(^V  )|  ^ 


2  -ArC/oHr-un-i 


WX  |  1 

1  0 


n(fi)  =  n(v) 
/»(//)  *  n(v) 


For  the  clutter  signal  we  note  again  the  one-to-one  correspondence  between  (p  and fp,  (p=xp  (fp), 
and  that  it  is  incoherent  over  <p. 

Following  an  argument  similar  to  that  preceding  (Cl 2),  where  we  showed  that  the  received  clutter 
power  spectrum 

Sc«  crc(fD)\F,(fD)?,  (E6) 

we  obtain  a  corresponding  received  clutter  voltage  vc  at  element  m,  time  delay  n 

~  >C (fo y*r(fn )eJfo )ej2xM'-(n~')T) ,  (E7) 


where  Oc  denotes  clutter  phase  angle. 

The  elements  of  the  clutter  matrix  Mc  are  thus 
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Kv  =  J  <?c(fD K(fo K (fn )* e-j2rf^-n^dfD  = 

=  J  o-c(/)  ko(/D)|2  e"yll(/t,)<r-0,,'r-"'»)  e-Wr>w*-w(ifD 


(E8) 


where  k (fD)  denotes  a  wave  vector  with  magnitude  2n  /  Xo  and  direction  corresponding  to/D. 
Finally,  adding  thermal  noise  PN  =  W  n0l  w"  to  the  clutter  gives  the  signal-to-clutter  +  noise  ratio 


pL 


— —  =  const.* 


pC  ,  pN 
out  out 


w'A  w” 
w'B  w 


(E9) 


where  matrix  A=M'  and  B  =  Mc +nQI . 

The  weight  vector  wmax  that  maximizes  this  ratio  of  quadratic  forms  is  obtained  by  using  a 
Singular  Value  Decomposition  (SVD)  to  compute  a  matrix  BU2  such  that 
Bu~ B]l~  =  B,  B~v 2  =  ( BV2y 1 .  The  desired  weight  vector  is  then  given  as 

wma  ,~(5-|/2emilx)*  (E10) 

where  enias  is  the  eigenvector  corresponding  to  the  largest  eigenvalue  of  (B~u2)+ A  B~U1  . 
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Appendix  F 

Satellite  Position  Accuracy 

The  accuracy  with  which  each  satellite  antenna  phase  center  is  known  determines  the  quality  of 
the  grating  lobe  pattern  formed  by  the  TechSat21  sparse  array.  Uncompensated  position  errors 
remove  power  from  the  array  grating  lobes  and  distribute  that  power  into  inter-grating  lobe  angles 
or  Doppler  frequencies.  Target  detectability  is  degraded  in  these  inter-grating  lobe  regions 
because  of  increased  clutter  power  at  Doppler  frequencies  competing  with  targets  at  these  same 
frequencies.  If  the  position  errors  become  large  enough,  performance  approaches  that  of  the 
random  sparse  array  discussed  in  the  last  section. 

We  assume  that  the  satellites  of  the  TechSat2l  array  are  in  a  single  orbital  plane  with  an  initial 
triangular  grid  described  in  Fig.  6b.  The  nominal  motion  of  each  satellite  is  governed  by  Kepler’s 
equations  for  a  spherical  earth  for  satellites  in  a  nearly  circular  orbit  (eccentricities  approximately 
10  s).  Deviations  from  the  nominal  satellite  positions  will  occur  due  to  other  forces  such  as 
residual  atmospheric  drag,  and  thruster  and  other  gravitational  forces.  We  assume  that  each 
satellite  has  onboard  thrusters  that  maintain  the  satellite  position  within  a  control  box  which  limits 
along-orbital-track  position  errors  to  be  less  than  ±e.  Here  e  is  of  the  order  of  meters.  Radial  or 
cross-track  control  errors  are  less  important  since  the  sparse  array  electronic  scan  control  can 
compensate  for  these  errors  at  a  constant  radar  range.  The  actual  satellite  position  with  respect  to 
other  satellites  in  the  array  is  measured  using  laser  ranging  to  provide  relative  positions  with  an 
error  significantly  less  than  a  wavelength  at  the  X  band  (10  GHz)  carrier  frequency  of  the  radar. 
Thus,  position  errors  of  the  satellites  are  modeled  as  deviations  from  the  nominal  positions  having 
a  random  but  known  component  on  the  order  of  meters  and  a  random  but  unknown  component  on 
the  order  of  1  O'3  meters  (=  0. 1  wavelength). 

The  random  but  known  position  errors  are  used  to  electronically  compensate  the  TechSat21 
sparse  array  to  form  a  beam  in  the  pointing  direction  of  the  array.  These  known  position  errors 
cause  grating  lobes  away  from  the  pointing  direction  to  degrade,  as  described  above.  This 
degradation  becomes  more  severe  for  grating  lobes  farther  away  from  the  pointing  direction  of  the 
array. 

To  assess  bounds  on  the  along-track  position  control  error  (±  e),  we  require  that  there  be  no  more 
than  a  3  dB  reduction  in  any  sparse  array  average  grating  lobe  within  the  main  beam  of  the 
satellite  pattern.  Clutter  power  from  increases  in  the  inter-grating  lobe  pattern  outside  the  satellite 
main  beam  will  be  reduced  to  a  two-way  sidelobe  level  of  the  satellite  pattern.  The  random  but 
known  position  error  of  each  satellite  is  assumed  to  be  uniformly  distributed  over  the  along  track 
control  box  from  -e  to  +e  and  independent  of  the  position  errors  of  the  other  satellites.  Using  this 
model,  the  reduction  in  average  grating  lobe  gain  of  the  one-way  sparse  array  pattern  due  to 
position  errors  was  computed.  Allowing  a  gain  reduction  of  3  dB  for  a  grating  lobe  at  the  first  null 
of  the  satellite  array  pattern  gives  the  following  bound  on  the  control  error  box: 

e  <  0.22L/B  . 

Here  L  is  the  along-track  satellite  aperture  dimension  and  B  is  the  azimuth  beam-broadening 
factor  [23]  associated  with  the  azimuth  aperture  distribution.  Here  B=1  corresponds  to  a 
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uniformly  illuminated,  maximum  gain  aperture  and  B=1.2,  1.36,  and  1.5  correspond  to 
Chebyshev  illuminations  with  -30,  -40,  and  -50  dB  one-way  sidelobes,  respectively.  Fig.  FI 
shows  this  result.  For  a  2m  by  2m  satellite  aperture,  the  control  box  varies  from  ±0.3  m  to  ±0.44 
m,  depending  on  the  beam-broadening  factor.  II  the  4  m"  satellite  array  aperture  is  conserved  to 
provide  adequate  signal-to-noise  but  distributed  as  a  4m  x  1  m  aperture,  the  decreased  azimuth 
beamwidth  increases  the  control  box  to  ±  0.6  m  or  ±  0.88  m.  This  reduces  the  satellite  thruster 
fuel  consumption  and  increases  mission  lifetime. 


Satellite  Alone  Track  Aperture  (m) 


Figure  FI.  Maximum  control  error,  c,  in  meters  vs.  the  satellite  along  track  aperture,  in  meters,  to  insure  that 
average  sparse  array  grating  lobes  within  the  aperture  main  beam  decrease  less  than  3  dB. 
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